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Axion cyclotron emissivity of magnetized white dwarfs and neutron stars 

M. Kachelriefi, C. Wilke, and G. Wunner 
Theoretische Physik I, Ruhr-Universitdt Bochum, D-44780 Bochum, Germany 

The energy loss rate of a magnetized electron gas emitting axions a due to the process e~ —>■ e~ + a 
is derived for arbitrary magnetic field strength B. Requiring that for a strongly magnetized neutron 
star the axion luminosity is smaller than the neutrino luminosity we obtain the bound Qae ^ 
for the axion electron coupling constant. This limit is considerably weaker than the bound derived 
earlier by Borisov and Grishina using the same method. Applying a similar argument to magnetic 
white dwarf stars results in the more stringent bound Qae ^ 9 • 
where T is the internal temperature of the white dwarf. 

PACS numbers: 14.80.Mz, 95.30.Gq, 97.60.Jd, 97.20.Rp 


I. INTRODUCTION 


The axion a is a pseudoscalar boson introduced by Pec- 
cei and Quinn (PQ) to solve the strong CP-problem in a 
natural way ||] . It results as a massless Goldstone boson 
from the spontaneous breaking of the PQ-symmetry by 
the VEV of a scalar fielT {( j )) = /pq. As noted first by 
Weinberg and Wilczek 0, the chiral anomaly of QCD 
induces, at low temperatures T < Aqcd ~ 200 MeV, 
a mass rria ~ ■^qcd//pQ the axion. The numerical 
value of the axion mass is given by 


0.60 X lO^GeV 
/pq 


eV. 


( 1 ) 


Besides the coupling to two photons through the chiral 
anomaly, the axion interacts with fermions by derivative 
couplings, which are all inversely proportional to the ax¬ 
ion mass. The interaction of axions with electrons is 
described by 


^ae = ( 2 ) 

where the coupling constant 

9ae — n Ce (3) 

/PQ 

depends on specific models through the effective PQ 
charge Cg. At tree level, the electron has Ce = 1/3 cos (3 in 
the DSFZ axion model, while it has Ce = 0 in the KSVZ 
model ||]. 

In invisible axion models, the axion mass rria is in prin¬ 
ciple arbitrary, but in fact severely constrained by astro- 
physical and cosmological considerations |^J^. The as- 
trophysical bounds on rria arise because axion emission is 
an additional energy loss mechanism for stars, and there¬ 
fore changes the stellar evolution. Since the axion cou¬ 
plings to matter are proportional to rria, axion emission 
is suppressed for small values of rria. Therefore astro¬ 
physics yields an upper bound for TOq. In ref. 1^, the 
bound rria ^ 9 • 10“^eV was derived for a DFSZ axion by 
studying the evolution of red giant stars. On the other 


hand, cosmology provides a lower bound for rria. Ax¬ 
ions would have been produced in the early Universe as 
coherent waves or by the decay of axion strings. There¬ 
fore they could constitute cold dark matter. Requiring 
that axions do not overdose the Universe yields the lower 
bound 10“®eV < TOq. 

It is the purpose of the present work to examine the 
emission of axions by electrons in the magnetic field 
5 of a strongly magnetized neutron star. The process 
considered, e“ ^ e“ -I- a, is very similar to cyclotron 
emission, e“ ^ e“ -I- 7, and, therefore, may be called 
axion cyclotron emission. This process was first ex¬ 
amined by Borisov and Grishina using a semiclassi- 
cal expression for the transition probability of the ele¬ 
mentary process e~ e~ + a, valid a priori only in 
the limit of low magnetic field strengths, B <C Bcr = 
TOg/e « 4.41 X 10^^ Gauss, and of ultrarelativistic elec¬ 
trons, E ^ rUg. By comparing the axion cyclotron emis¬ 
sion with neutrino cyclotron emission e~ e~ -\- v + v, 
they concluded that the axion electron coupling constant 
is bounded from above by Qae % 5 • 10“^'*. If this re¬ 
sult, corresponding to rria ^ 6 • 10“^eV for Ce = 1, 
could be confirmed or improved, the window in the ax¬ 
ion mass range would become very narrow, or even be 
closed. Consequently, it is imperative to recalculate the 
axion cyclotron emissivity and to extend the results ob¬ 
tained in ref. Q to arbitrary magnetic field strengths and 
electron energies. We derive the probability of the ele¬ 
mentary process e“ ^ e“ -I- o taking the magnetic field 
exactly into account, and then numerically evaluate an 
expression for the axion emissivity which is exact in the 
limit of a completely degenerate Fermi-Dirac gas. As the 
main result of this work we obtain the axion cyclotron 
luminosity of a degenerate star with arbitrary core mag¬ 
netic field strength B. Although we confirm the analyt¬ 
ical approximations of ref. Q| in the limit B <?; Bcr and 
E ^ rUg, the bound derived by us for the axion electron 
coupling constant applying the energy-loss argument to 
neutron stars is three order of magnitudes weaker. The 
main reason for this is that all conventional neutrino 
emission mechanisms such as the modified URCA pro¬ 
cess, which we take into account, were neglected in ref. 
iQ. Finally, we consider the energy loss of white dwarf 
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stars due to axion cyclotron emissivity. Requiring that 
for a white dwarf star the axion luminosity is smaller 
than the photon surface luminosity we obtain the bound 

II. ELEMENTARY PROCESS 


The time integration and two of the momentum inte¬ 
grations yield the usual conservation laws, 


E = E' + Lu 

(15) 

Py=Py+ ky 

(16) 

+ 

II 

h} 

(17) 


We start by calculating the probability that an electron 
with polarization r = ±1 occupying an excited Landau 
level N emits an axion a with momentum k = {uj,k) and 
undergoes a transition to the Landau level N'. (Quan¬ 
tities describing the final electron are primed; compare 
also the similar calculations of e“ —> e“ -I -7 in ref. j9| for 
more technical details.) The ^'-matrix element for this 
process, 


%= , ( 4 ) 

follows directly from the Lagrangian (I) if we assume 
that there exists only one Goldstone boson. 

We choose the homogeneous magnetic field in z direc¬ 
tion, B = Bcz, and as wave functions for the electron 


while the x integration results, up to a phase factor, in 
Laguerre functions I, 


In',n{k) 


,JN-N')/2 

V N\ 



(18) 


for N > N'. The argument of the /-functions is given 
by K = (fcsin0)^/(2ei?), 9 is the angle between B and 
fc, and the Laguerre polynomials L" (x) are defined as in 
ref. [^. Using Eq. dil) and we obtain the energy 
Lo of the axion as a function of 9 

uj{9) = — pz cos 9 — — pz cos 9)^ 

-2{N - N')eB 9 , (19) 


tpix) 


(eR)^/^ g-i(£^t-Pyy-p^2) 

yj LyLz 2V EEq 


( Ai \ 

A 2 1 

Ki (j)N-l{0 I 

V K2 4>n{0 / 


( 5 ) 


Here, ipNiC ~ ^ ^^e normalized eigenfunctions 

of the harmonic oscillator, and we have introduced the 
abbreviations 


where, as in the following, we have neglected the axion 
mass 

The decay width of an electron in the Landau level N 
into an Landau level N\ 

{k) = p Y, \S±fT^d^k, ( 20 ) 

P'y X ,r' 


Ai = 7/(1 + r) + A (1 - r) = nt 
A 2 = A (1-h r)-f i/(1 - r) = n; 
ATi = Z (1 -t- r) - r (1 - t) = Ml* 

7/2 = r (1 -f r) - Z (1 - r) = M* 


and 


( 6 ) 

( 7 ) 

( 8 ) 
(9) 


where ± denotes the (in general different) decay widths 
of Landau states with polarization r = ±, becomes 



wsin 0 \M\^ 

162 EEoE'q{E' -p'^ cos 9) 


( 21 ) 


7/ = y/(E + Eo)(Eo + Trie) 


A = -ip. 


Eq - me 
E Eq 


Z =Pz 


Eq + rUe 
E + Eq 


r = i^/{E + Eq){Eq - me) 

Eq = v'E2 - pI = ^ml + 2NeB . 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


with 


M = (n'lT/i - M(Ai)/jv'-i.iv-i(«) 

-(n^7/2 - M'A2)/Ar-,jv(Aj) 


+ + 


and aae = 


( 22 ) 


^We will see below that rria <C ta < E_f, where Ep = 0{me) 
is the Fermi energy of the degenerate electron gas. 
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III. AXION EMISSIVITY 

The luminosity La emitted by M electrons occupying 
the volume V due to the process e~ e~ + a is given by 


f 


dpz g{E)S « g^Ep) 

. Pz 


-I' 


(30) 




T 


(23) 


A.A' 


where the summation index A indicates the set of quan¬ 
tum numbers A = {N,T,py,pz}, 

S = f{E) [1 - f{E')] (24) 

and f{E) are Fermi-Dirac distributions functions 

fiE) = [exp{(3iE-p)) + l]-\ (25) 

For the electron density rig = -^/{ExLyL^) and temper¬ 
ature T = the chemical potential p is fixed through 


and the Fermi energy E-p can be approximated by Ep « 
p. Now the summation over the discrete Landau levels 
N breaks up at A„iax = int [(i?| — m^)/{2eB)], while 
the momentum p^ of the initial electron is determined 
by pz = \JE'^ — Eq. When {E^ — m1)/i2eB) is an in¬ 
teger, the axion emissivity diverges because Pz is zero or 
because, more physically speaking, the density of initial 
states is infinite. This can be remedied by taking into 
account the finite life time of excited Landau states (cf. 
ref. 1^) but is unimportant for our purposes. 

Although the two infinite sums over N and N' have 
been replaced by finite sums, for low magnetic field 
strengths or high densities it is still necessary to com¬ 
pute Laguerre polynomials with high index. This can 
be avoided in the semiclassical, ultrarelativistic case 
{B/Bcr I, E me) by the use of a Bessel function 
approximation for the /-functions H) 




^ _ pOC t 


(26) 

Here the factor eB L^Ly/{2 tt) takes into account the Py 
degeneracy of the Landau states. For a nearly degener¬ 
ate Fermi-Dirac gas, the Pz integral in Eq. (|2m can be 
approximated (in a similar manner as in ref. |ll|) by 


with 


and 


Ko= [VN- VW 


Ei ^2 


6 


p^ - El 


O (T^) 


2 

^ ( 2 
z = 7;{ko 


NN') 


1/4 


1 


KoJ 


3/2 


(31) 


(32) 


(33) 


(27) 


This asymptotic expansion is only valid for k \ kq, i.e. 
0 « 7r/2 is a necessary condition for its applicability. 
Therefore we evaluate dT for Pz = 0 and then carry out 
a Lorentz transformation to Pz = \JE^ — El , viz. 


where Amax = int [(/i^ - ml)J{2eB)]. 

The axion cyclotron emissivity Sa of electrons, i.e. 
the rate of energy loss per volume due to the process 
e“ ^ e“ + a, follows from Eq. (p^) as 


ea = 


eB 


E E E 

N=1 N’<N T=±' 


dpz / dO uj- 

> Jo 


dT 


N^N' 


{k) 


dO 


S. 


COS0Q = 


dTlf) = dTl 


E cos 0 — Pz 
E — Pz cos 9 


Eq 

E — Pz cos 9 


(28) 


dT{k) = ^ dFo(fco). 


(34) 


(35) 


(36) 


The numerical evaluation of the terms in Eq. (28) be¬ 
comes cumbersome already for moderate N. Therefore 
we take advantage of the degeneracy of the electron gas 
inside a white dwarf or neutron star and employ an ap¬ 
proximation commonly used in calculations of neutrino 
emission rates: With the identity 


Here quantities evaluated in the rest frame of the initial 
electron are denoted by the subscript 0. Insertiiig the ap - 
proximation for S into Eq. (|2|) and using eq. (P^-(p^, 
we obtain as final result for the axion emissivity in the 
limit of a degenerate electron gas 


5 = 


1 


= [/(^) - 


(29) 


cB V -V V -V V -v 
“ ^ 9 ^ 2^ 2^ 2^ 




2tt^ 


and f{E) « 9{E — p), f{E') « 9{E' — p), the integration 
of a slowly varying function g{E) of E over pz becomes 


N=1 N'<N T=±-^0 

Eo 


d9 sin 9 — 

Pz 


E — Pz cos 9 


2 ^n^n' 

di j_Q 


(fc) 


UI 


dOo 


e/3<^ - 1 ■ 


(37) 
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Finally we examine the classical limit p], <C m\ and 
B Bcr- Inserting the approximations E « irie + 
NeB/irie, E' « me + N'eB/me +p'z^l (2me), and uj{6) ft! 
(TV — N')eB/me in \M\'^, we see that \M\^ is of order 
0{B). Moreover, since — 1]“^ ft T/uj, the axion 
emissivity is proportional to B'^T in this limit. If the 
condition N^B/Bcv ^ 1 is also valid, we can approxi¬ 
mate the /-functions by 






^/W\ 


(38) 


IV. NUMERICAL RESULTS 


treatment of a magnetized Fermi-Dirac gas. For rela¬ 
tively low Neaax, the emissivity decreases exponentially 
until for 2eB > p? — all electrons populate the Lan¬ 
dau ground state (for T = 0 ) and the emissivity becomes 
zero. 

To obtain a bound for the coupling constant pae, we 
now compare the axion luminosity with other processes 
important for the cooling of neutron stars. A process spe¬ 
cific to magnetized neutron stars is the cyclotron emission 
of neutrino pairs by electrons, e~ e~ + v+D. Kaminker 
et al. Q obtained 

= (9 X 10“erg cm-’ s-' ) (Ajj) 

(40) 


A. Neutron stars 

First, we will establish that Eq. (|37| ) derived by us for 
general B and Ue has the correct semiclassical limit for 
E^ me and B ^ B^ as given in ref. 0, 




cl 

a 


= 8.67 • IQ-^glem 
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(39) 


Since for a relativistic degenerate Fermi gas and 
1/3 /Tp 2eB the usual relation pp = -^STr^ne can 
be used 0], this approximation is applicable for rig 2> 
mg/(37r^) ft 5.84 • 10^®cm“^. In Fig. we com¬ 
pare our numerically accurate results with the approx¬ 
imate results for = ea/aiev as functions of B/Ba 
for He = lO^^cm”^, and for temperatures T = 0.04me 
and T = O.lmg, respectively. (We use the emissivity 
scaled by aiev = Q!oe/(5.75 • 10“^^) to obtain a quantity 
independent of TOo-) One recognizes that the exact re¬ 
sults (crosses and diamonds) reproduce the correct semi¬ 
classical B and T dependence (straight lines) already for 
B/Bct ^ 1- To demonstrate the saw-tooth-like behavior 
caused by the root singularities in the number density 
of the Landau states, an amplification of a small excerpt 
of Fig. 0 is displayed in Fig. ||. Note that these singu¬ 
larities are cancelled out on average in the semiclassical 
approximation. 

Fig. I shows the magnetic field dependence of the 
axion emissivity Sa for ne = 8.4 • 10^®cm“^ and T = 
0.04me ft 2.4 • 10®K. Assuming an ideal mixture of nu¬ 
cleons and electrons in //-equilibrium, this value of the 
electron density rig corresponds to the nuclear density 
po = 2.8 • . It can be seen that the semiclas¬ 

sical approximation produces reasonable results even for 
B/Bct ~ 20, and not only for B/Bcr ^1. A more 
suitable criterion for the applicability of Eq. (^) is 
A^max 1, or, equivalently, for an ultrarelativistic de¬ 
generate electron gas, Ep ^ 2eB. Except for a fac¬ 
tor of three, this criterion is the same as the condition 
derived in ref. for the validity of the semiclassical 


for the energy loss rate due to this process. The bound 
Pae ^ 5 • 10 “^"^ of ref. was derived by requiring 
^ Su.cyc under conditions appropriate for the outer 
crust of a neutron star. Since only the total luminosity of 
stars is observable, a more suitable - but weaker - crite¬ 
rion is to demand that the axion luminosity Ca is smaller 
than the total neutrino luminosity L^,. Other energy loss 
processes contributing to Ci, are neutrino bremsstrahlung 
in the Coulomb field of ions iMl]’ 

Eiong = (2.1 X 10"°erg cm'^ g-^)^ ^ ( 1^) 
and the modified URCA process, 


£urca = (2.7 X lO^^erg cm'^ g-^) ( ^ ) ^ 

m — 


TOg 




T 


lO^Kj 


(42) 


The value jA is in the range 1 — 10 [^, and TO*,m* 
are effective masses of the Landau theory. In the fol¬ 
lowing we use jA = 5 and Wn / m* = mp/m* = 1. 
Recently, Pethick and Thorsson Q recalculated the neu¬ 
trino bremsstrahlung process. Instead of treating the 
interaction of the electrons with the ions immersed in 
a lattice in first-order perturbation theory, they used 
Bloch wave functions for the electrons and found that 
bremsstrahlung is exponentially suppressed at low tem¬ 
peratures, T < mg. 

As a model neutron star we choose a neutron star 
with radius R — 10 km, constant density p equal to nu¬ 
clear density po and temperature T = 0.04mg. Since 
the neutron star observations by ROSAT indicate that 
the standard cooling scenario without exotic matter de¬ 
scribes the cooling curves , we assume the absence of 
a pion condensate. Young pulsars have magnetic dipole 
fields lO^^G < Sdipoie ^ 3 • lO^^G. Since the magnetic 
field strength inside the neutron star is generally assumed 
to be stronger than Rdipoie (cf. e.g. | 20 |), we choose 
B = 0{1- lORcr). 
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First we consider the case where the core of the neutron 
star is not superfluid. Then the total neutrino luminosity 
is Lv = R^ision + eurca) « 1-3 - 8.9 • 10^®erg/s. 
Here the first number takes into account the exponen¬ 
tially suppression found in ref. ||^ by a factor 0.02 for 
Eion, while the second number uses Sion given by eq. (^). 
The luminosity due to cyclotron emission of v pairs (eq. 
(^)) could be neglected in both cases. For T = O.OdrUe, 
we read off the axion emissivity £a ~ lO^^erg cm“^s“^ 
from Fig. ^ and obtain for the axion luminosity Ca = 
4 • 10^® aiev erg/s. Requiring £a < we arrive at the 
bound gae ^ 5 • 10“^^ — 1 • 10“^°. 

Next we comment on the case that in the core of the 
neutron star superfluid protons form a type II supercon¬ 
ductor 1^] . Then the initially homogenous magnetic field 
B becomes confined to quantized flux tubes (Abrikosov 
fluxoids). Each fluxoid carries an elementary magnetic 
flux quantum (pQ = -nje and has the field profile 


where A is the penetration depth of the magnetic field 
and ATo is a modified Bessel function. The above for¬ 
mula for B{r) is valid if the proton-proton correlation 
length is much smaller than the distance r from the flux¬ 
oid axis. The fluxoids are separated by the mean dis¬ 
tance dp = [2(/)o/(-\/3H)]^/^. We can approximate the 
mean axion emissivity (Eq) by the average over the lat¬ 
tice of fluxoids with a quasi-homogenous magnetic field 
given by eq. (|3[), 


2 

{Sa) ~ J drr£a{r), (44) 

if the coherent interaction length Zjnt is much smaller than 
the characteristic length scale A on which the magnetic 
field varies. Typical values for A are A « 30 fm 
while Zint is given by the uncertainty principle, hut,x,y ^ 
Xj^feB and by lint,z ^ I/Pz- Thus Zint has the same order 
of magnitude as A and the quasi-uniform approximation 
is not well justified. Therefore the approximation ( ^ ) 
should be seen only as a crude estimate. For initial mag¬ 
netic field strengths B = Bcr and B = lOHcr we obtain 
(Eq) ~ 3-lO^^erg cm“^s“^ and (Eq) « 3-10^®erg cm“^s“^, 
respectively. This means that axion cyclotron emissivity 
is only weakly suppressed by the superfluidity of pro¬ 
tons, while the neutrino emissivities & and (^) are 
suppressed by a factor of exp(—2/3A„), where A„ is the 
energy gap at T = 0. On the other hand, new energy 
loss mechanisms such as neutrino emission by scatter¬ 
ing of electrons on fluxoids could become important . 
Because of the uncertainties in both the axion and the 
neutrino emissivities, we do not consider it appropriate 
to try to derive a rigorous bound for gae in the superfluid 
case. 


B. White dwarfs 


We now apply the same line of arguments to magnetic 
white dwarfs. The internal temperature T of a white 
dwarf with mass M is related to its surface photon lumi¬ 
nosity by §J2|] 


M 


3.3- 10“^ 


erg 

gs 



(45) 


About 5% of white dwarfs have surface magnetic fields 
in the range 10® — 10®G. Although the internal magnetic 
fields strengths are assumed to be in the range 10® — 10^® 
G, they could approach even 10^^ G fl^ . 

The mean density of a white dwarf depends on its 
total mass and chemical composition but nevertheless 
a typical mean density is 10®g/cm® corresponding to 
Ue ~ i ■ 10^®cm“®. Unfortunately, for this value of Ue 
none of the two approximations for the /-functions is 
applicable. Therefore, we compute £a for densities low 
enough such that the use of the approximation Eq. (|^ ) 
is possible. In Fig. Eq is shown for rig = 10^®cm“® and 
the temperatures T = 10®K, T = lO’^K and T = 10®K, 
respectively, together with the fit function 


;rfit 


= 2.6 


erg 

cm®s 



(46) 


The agreement with the assumed B^T behavior is good. 
We want to give a simple argument in order to decide 
if eq. (M) over- or underestimates the true emissivity 


because of the use of the low density Ue = 10^®cm 


- 3 . 


The ultrarelativistic approximation eq. 


yields, for 


the values B = 10®G, T = lO^K and rie = 3 • 10^®cm“®, 
a much higher emissivity tha n Therefore the non- 
relativistic approximation eq. (46) can be used as a lower 
bound for the true emissivity, if one assumes a smooth 
transition between the non-relativistic and the ultrarela¬ 
tivistic asymptotic behavior. 

We now require that magnetized white dwarfs do not 
emit more energy in axions than in photons, i.e. that 
£j/M < Ea/p, and obtain 


flae < 9 • 10 


-13 


/ j. \ Q y 

Vio^k; uQioGy 


(47) 


Since this bound is quite sensitive to B and the knowl¬ 
edge of the internal magnetic field strengths of white 
dwarfs is poor, it is hard to derive a precise bound for 
gae- However, since Eq oc T, compared to C-y oc ax¬ 

ion emission becomes more and more important during 
the cooling history of magnetized white dwarfs. There¬ 
fore the fraction of magnetized white dwarfs among all 
white dwarfs should diminish drastically for low enough 
luminosities. In Table I, we compare the fraction of 
strongly magnetized white dwarfs with the fraction of all 
white dwarfs in three different temperature bins [p|,p4[. 
The distribution of hot strongly magnetized white dwarfs 
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adapted from Ref. [p4[ spans approximately three order 
of magnitudes in surface dipole field strength, Bdipoie = 
3 • 10® —10® G, and has a maximum at Bdipoie ~ 3 • 10^ G. 
It consists of only 18 stars, so some caution in interpret¬ 
ing the data is appropriate. Nevertheless, the fraction 
of magnetized white dwarfs in the last temperature bin 
is considerably diminished compared to the total white 
dwarf population. One possible explanation for this could 
be the additional energy loss of magnetized white dwarfs 
due to axion cyclotron emission. 


V. SUMMARY 


We have derived the axion emissivity of a magnetized 
electron gas due to the process e~ —^e~+a for arbitrary 
magnetic field strengths B. Gomparing the exact emis¬ 
sivity with the semiclassical approximation obtained ear¬ 
lier, we have shown that this approximation is not only 
valid for B/ 3^ 1, but also under the more general 

condition Ep ^ 2eB. Applying axion cyclotron emission 
as an additional cooling mechanism to neutron stars and 
requiring that the axion luminosity is smaller than the 
neutrino luminosity we could constrain the axion electron 
coupling constant to a value as small as pae ^ 0(10“^®). 
This bound is three orders of magnitude weaker than the 
bound derived in ref. 0 considering the same process. 
This discrepancy is not caused by differences in the cal¬ 
culation of the axion emissivity but by different assump¬ 
tions about neutron star cooling. In the case of white 
dwarfs we could derive the more stringent limit eq. (47). 
We have noted that the lack of low-temperature magne¬ 
tized white dwarfs could be interpreted as signature of an 
additional energy loss due to axion cyclotron emission. 


ACKNOWLEDGMENTS 


During the initial phase of this work, MK was sup¬ 
ported by grants from Deutscher Akademischer Aus- 
tauschdienst and Land Nordrhein-Westfalen. 


[4] E. W. Kolb and M. S. Turner, The Early Universe 
(Addison-Wesley: Redwood City 1990). 

[5] G. Raffelt, Stars as Laboratories for Fundamental 
Physics (University of Chicago Press: Chicago 1996). 

[6] G. Raffelt and A. Weiss, Phys. Rev. D51, 1495 (1995). 

[7] A. V. Borisov and V. Yu. Grishina, JETP 79, 837 (1994). 

[8] A. A. Sokolov and I. M. Ternov, Synchrotron Radiation 
(Akademie-Verlag: Berlin 1968). 

(91 H. Herold, H. Ruder, and G. Wunner, Astron. Astrophys. 
115, 90 (1982). 

[10] M. Abramowitz and I. A. Stegun (eds.). Handbook of 
Mathematical Functions (Dover: New York 1970). 

[11] L. D. Landau and E. M. Lifshitz, Statistical Physics, p. 
162f (Pergamon Press: London 1959). 

[12] V. Ganuto and H.-Y. Chiu, Phys. Rev. 173, 1220 (1968). 

[13] A. D. Kaminker, K. P. Levenfish, and D. G. Yakovlev, 
Sov. Astron. Lett. 17, 450 (1991). 

[14] B. L. Friman and O. V. Maxwell, Astrophys. J. 232, 541 
(1979). 

[15] S. L. Shapiro and S. A. Teukolsky, Black Holes, White 
Dwarfs, and Neutron Stars (Wiley: New York 1983). 

[16] N. Straumann, General Relativity and Relativistic Astro¬ 
physics (Springer-Verlag: Berlin 1991). 

[17] J. W. Negele and D. Vautherin, Nucl. Phys. A207, 298 
(1973). 

[18] C. J. Pethick and V. Thorsson, Phys. Rev. Lett. 72, 1964 
(1994). 

[19] W. Becker, Ann. N.Y. Acad. Sci. 759, 250 (1995). 

[20] C. Thompson and R. C. Duncan, Astrophys. J. 408, 194 
(1993). 

[21] For a review see e.g. J. A. Sauls in: H. Ogelman and E. 
P. J. van den Heuvel (eds), Timing Neutron Stars, p. 457 
(Kluwer Academic Publishers, Dordrecht 1989). 

[22] A. D. Kaminker, D . G. Yakovlev and P. Haensel, preprint 
astro-ph/970215^ , to appear in Astron. Astrophys. 1997. 

[23] L. Mestel, Mon. Not. R. Astron. Soc. 112, 583 (1952). 

[24] G. D. Schmidt and P. S. Smith, Astrophys. J. 448, 305 
(1995). 


[1] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 
(1977); Phys. Rev. D16, 1791 (1977). 

[2] S. Weinberg, Phys. Rev. Lett. 40, 223 (1978); F. Wilczek, 
Phys. Rev. Lett. 40, 271 (1978). 

[3] M. Dine, W. Fischler, and M. Srednicki, Phys. Lett. 
B104, 199 (1981); A. P. Zhitnisky, Sov. J. Nucl. Phys. 
31,260 (1980). J.-E. Kim, Phys. Rev. Lett. 43, 103 
(1979); M. A. Shifman, A. I. Vainshtein, and V. I. Za¬ 
kharov, Nucl. Phys. B166, 260 (1980). For a review 
of different axion models and references see J.-E. Kim, 
Phys. Rep. 150, 1 (1987). 


6 







TeffllO^K] Fraction of all WDs Fraction of magnetized WDs 


40-80 

1% 

0% 

20-40 

23% 

50% 

12-20 

76% 

50% 


TABLE I. Fraction of all white dwarfs (WDs) and of mag¬ 
netized WDs in three temperature bins. 
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FIG. 1. Axion emissivity ia due to e —> e + a as a function of BfBa for Ue = 10^^ cm ^ and for T = 0.04me (bottom) 
and T = O.lme (top); lines: asymptotic formula, crosses (+) and diamonds (o): exact results). 


FIG. 2. Amplification of a small part of Fig. showing the saw-tooth-like behavior of the exact axion emissivity (for 
T = O.OdrUe). 
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FIG. 3. Axion emissivity Sa due to e —> e + a as a function of B jBcv for rie = 8.4 • 10^®cm ® and T = 0.04me (dashed 
line: asymptotic formula, diamonds: exact results). 


FIG. 4. Axion emissivity ia due to e —> e + a as function of B/B^ for Ue = 10^®cm ® and T = 10®K, T = lO’^K, 
T = 10®K (lines: fit function Eq. (44), points: exact results). 
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Figure 1: Axion emissivity Sa due to e“ ^ e“ + a as a function of BfBcr for 
Tie = and for T = 0.04me (bottom) and T = O.lme (top); lines: 

asymptotic formula, crosses (+) and diamonds (o): exact results. 
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Figure 2: Amplification of a small part of Fig. 1 showing the saw-tooth-like 
behavior of the exact axion emissivity (for T = 0.04me). 
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Figure 3: Axion emissivity Sa due to e“ ^ e“ + a as a function of B/Bcr 
for He = 8.4 • and T = 0.04me (dashed line: asymptotic formula, 

diamonds: exact results). 
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Figure 4: Axion emissivity due to + a as function of B/Ba for 

rig = and T = 10®K, T = lO^K, T = 10®K (lines: fit function Eq. 

(44), points: exact results). 
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